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ABSTRACT 


Kirchhoff 's formula for radiation from a closed surface has been used 
recently for prediction of the noise of high speed rotors and propellers. 
Because the closed surface on which the boundary data are prescribed in these 
cases is in motion, an extension of Kirchhoff's formula to this condition is 
required. In this paper such a formula, obtained originally by Morgans for the 
interior problem, is derived for regions exterior to surfaces moving at speeds 
below the wave propagation speed by malcinq use of some results of generalized 
function theory. It is shown that the usual Kirchhoff formula is a special 
case of the main result of the paper. The general result applies to a deform- 
able surface. However, the special form it assumes for a rigid surface in 
motion is also noted. Some possible areas of application of the formula to 
problems of current interest in aeroacous tics are discussed. 
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INTRODUCTION 


Kirchhoff's formula, published in 1882, is used in the theory of diffraction 
of liqht and in other electromagnetic problems [1-2] • It also has many appli- 
cations to problems of wave propaaation in acoustics [ 3 ] • One of the novel 
uses of this formula was proposed by Hawkinqs for predicting the noise of high 
speed propellers and helicopter rotors [4]. His idea involves surrounding the 
rotating blades by a closed surface S which moves with the forward speed of the 
machine. Inside this surface, nonlinear aerodynamic calculations are carried 
out which give the blade loads and the pressure history and its spatial and 
temporal derivatives on the surrounding surface S. In the exterior of this 
surface, it is proposed that a formula similar to Kirchhoff 's be used to calcu- 
late propagation of sound in terms of these surface values. This necessitates 
extension of the Kirchhoff formula to apply to moving surfaces. In 1930, such 
an extension was derived for the interior problem by Morgans [5]. His analysis 
was lengthy and somewhat complicated, and, in addition, it does not seem to be 
well known among acousticians. For these reasons, a modern derivation of a 
generalized Kirchhoff formula applicable to moving surfaces is presented 
here. It is very likely that specialized versions of this formula have been 
rederived by many researchers in the context of specific applications. One 
such case known to the authors is discussed in the last section of the current 
paper. 

In the following section the method used in the derivation of the main 
result is illustrated by deriving the Kirchhoff formula for radiation into the 
reuion exterior to a stationary surface S. The method involves application of 
the wave operator to a function $ which is set equal to zero outside the domain 
of interest. The function $ is thus discontinuous across the surface S. The 
derivatives in the wave operator are considered as generalized derivatives 
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[6-9] which are defined everywhere* The resulting expression then contains 
terms which involve the jumps in $ and its derivatives multiplied by Dirac 
delta functions whose supports are on S. These terms can be viewed as source 
terms for a wave equation which is valid in the entire unbounded space. 
Kirchhoff's formula is then derived using the Green's function of the wave 
equation for unbounded space. 

In the main section of this paper, the above method is extended to a 
smooth movinq surface* It will be observed that essentially no major difficul- 
ties arise in this case. One thus avoids the conventional approach which can 
become complicated in the case of moving surfaces. The resulting general 
Kirchhoff formula is then written for piecewise smooth rigid surfaces. It is 
also shown that the usual formula for a stationary surface is a special case of 
the qeneral result presented here. 

In this paper the wave propagation problem is treated as an acoustic prob- 
lem and the terminology of acoustics is used. The main result is of course 
applicable to many other wave propagation problems. Immediate applications to 
the problem of the prediction of the noise of high speed propellers and rotors 
are expected. 


EXAMPLE OF THE MATHEMATICAL APPROACH 
In this section the mathematical method of deriving the main result is 
illustrated by utilizing it to obtain the Kirchhoff formula for a stationary 
surface [1-3]. Consider the closed and bounded smooth surface S described by 
f(x)=0 such that f > 0 in the exterior of this surface. Assume that # and its 
first derivatives are continuous in f > 0 and that $ satisfies the wave 
equation 


i lit . V 2* =r-p* , o ( i ) 

C 2 at 2 
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in the exterior of surface S r and that $, 4» and <J> are qiven on S for 

t n 

te(-®,»). Extend $ to the interior of S(f < 0) by assuming that $ = 0 inside 
S. To reduce confusion between $ and this extended function, the latter 
function <fr(x,t) is defined as follows: 

f > 0 (exterior of S) 
f < 0 (interior of S) 

The function $ is a discontinuous function whose derivatives do not exist 
on f =0 • One can, however, define generalized derivatives for such a function 
[6-9]. In this paper operators involving generalized derivatives will be de- 
noted by a bar over the operator symbol . It follows from the definition of 
<1> that wherever ordinary derivatives of ? exist (i.e., interior or exterior to 
S ) one has 


*(x,t) = 


f $(x, t ) 

1 cT 


□ 2 * = 0 . 


(3) 


However, one can show, as is done below, that 


JL !_! - 7 2 * * 0 

c 2 3t 2 


(4) 


It will be seen that what appears on the right side of Eg. (4) leads to the 
Kirchhoff formula. For the reader who is not familiar with generalized func- 
tion theory, it is mentioned that one requires considerable conceptual develop- 
ment in goinq from Eq . (3) to Ea. (4). Examples of applications given in 
references [7] and [8] mav help in understanding the following steps. 
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One can always define the surface f=0 in such a way that jvf|=1 on this 
surface. This is done by taking df=dn where dn is distance from the surface 
f =0 alonn its local normal. With this definition one has 7f = n where n is the 

~ /V 

local unit normal pointing into the exterior region. Then, using generalized 
derivatives 18 ], one can write at once 


at 


9 $ 

at 


-- + 


3t 


6(f) 



(5-a) 


111 = 111 

at 2 at 2 


( 5-b) 


7# = V* + * n 6(f) , 


(5-c) 


? 2 $ = v 2 * + $ 5 (f) + 7 • [ $n 6 ( f ) ) , 

n ~ 


(5-d) 


where in Eq . (5-c) the function <J> on the right side is defined as 


$ = lim $(x,t) ( 6 ) 

f+0 + 


The function $ in Eq . (5-d) is the normal derivative of and is defined 
n 

similarly. The Dirac delta function is denoted 5(f)* From Eq. (5) one obtains 


□ 2 ? =n 2 t - * 6(f) - 7»[*n6(f)] = -* 6(f) - 7»(*n6(f)] • (7) 

‘ ' n ~ n ~ 


Here, Eq . (3) is used in the last step. This completes the determination of 
the right side of Eq . (4). 

Equation (7) is valid in the entire unbounded space. One can now use the 
Green’s Function for the wave equation in unbounded space, 5(g)/4Trr, where 
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q=T-t+r/c and r= | x— y | . The space-time variables of the observer and source 
are denoted (x,t) and (y,t)f respectively. The speed of sound is c. Equation 
( 7 ) then yields 

$ $n 

4w*(x,t) = - / — 6(f)fi(g)dydx - V •/ — 6(f )6(g)dydx. (8) 

The space inteqrals in Eq. (8) are over the entire unbounded space and the time 
inteqral is over (-»,t]. The divergence operator on the right acts on the 
variable x. r Jse the fact that the volume element d£ = dfdS, where dS is the 
element of area of the surface S: f = 0. Then let t + g noting that the 

Jacobian of the transformation | 3 1 / 3g | = | 3^ / 3 T | 1 is unity. Integration over 

f, by virtue of the delta function 6 (f), restricts the expressions to f = 0. 
Subsequent integration over q restricts the expressions to g = 0, which intro- 
duces the retarded time. Equation (8) thus becomes 

4ir*(x,t) = - / — — dS - 7 •/ n dS (9) 

^ j r J r ~ 

f =0 f =0 

Here the square bracket stands for evaluation of the functions at the retarded 
time; i.e., [$>] = #(y,t-r/c) . 

Taking the gradient operator inside the second integral of Bq. (9) and 
collecting terms of 0(1 /it ) and 0(1/r^), one gets 

4tt$(x, t) f > 0 

4tt$(x , t ) = { 

0 f < 0 


- / 

f=0 


[c 1 $ cos 9 - # ] 
T n 


dS + 


/ 

f=0 


[ $]cos9 

r2 


dS 


( 10 ) 
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Here <t>^ = 94 >(£,t)/9t and 0 is the anqle between n (pointinq into the exterior 
region) and the radiation direction r=x-y. For the exterior region f > 0, this 
is Kirchhoff's formula for a stationary surface S. The present method auto- 
matically shows that for the observer inside the surface S the inteqrals in 
Kirchhoff's formula yield a zero (null) field. This means that and 

are not independent on the surface S. Equation (10) is an identity which may 
be utilized in various ways. As it stands it may be considered as an integral 
representation of 4> at points exterior to S in terms of surface quantities. 
Another interpretation follows if the observer point x is taken on S itself; in 
that case, Eq . (10) becomes an inteqral equation governing $ on S. This latter 
interpretation is the basis of modern boundary element methods. It is seen 
that the integrals of Eq . (10) remain valid for a piecewise smooth surface S. 

It is noted that in the usual applications of Eq. (10) $ is taken as a 
complex harmonic function of time so that the temporal part of <J> disappears in 
Ea. (10); the resultinq Kirchhoff formula pertains to the complex amplitude of 
$ [1-3], The amplitude satisfies Helmholtz' equation. In this case one could 
begin the above procedure from a Helmholtz equation rather than from the wave 
equation as is done above. 


THE GENERAL KIRCHHOFF FORMULA 

In this section, the Kirchhoff formula for moving surfaces will be derived. 
Let f{x,t) = 0 describe the moving surface S which is assumed to be piecewise 
smooth. The function f is defined as before such that f > 0 outside S and Vf 
= n on f = 0. The function $ and its first derivatives are assumed to be 
continuous and $ satisfies Ea . (1). This function is again extended to $ by 
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Eq. (2) so that Bq . (3) is satisfied both inside and outside S. The terms on 
the right side of Eq. (4) are found by the same technique as presented above. 
Taking the qeneralized derivative of $ with respect to time qives 


— 'V 

3* 

3t 


3$ 

3t 


= — + * 


£ «*> 


3 * 

at 


♦ v 6(f) 
n 


( 11 ) 


where v n = -3f/3t is the local normal velocity of S with respect to the undis- 
turbed medium. Dif ferentiatinq both sides of Eq. (11) once more yields 


HI = 32 ? 

at 2 at 2 


6(f) 

t (x ) n 


3t 


[*v 6(f)] 
n 


( 12 ) 


where = 3<t/3t|^ . The qeneralized Laplacian of $ is qiven by Eq. (5-d) as 

before. From this and Eq . (12), one obtains 


□ 2 ? =-($ + ! M $ . ) 6( f ) - -4- [M $6( f ) ] - V* [ ♦ n 6(f)] 

ncnt(x; c3tn ~ 


(13) 


where M n = v n / c I s tlie local normal Mach number on S, In the following, 
partial derivatives with respect to time holdinq different sets of space vari- 
ables fixed will appear* The notation introduced above (e*g. will be 

maintained wherever necessary in order to prevent confusion. The source terms 
in Eq. (13) can be written in other forms by performing the temporal and spa- 
tial derivatives in the second and the third terms, respectively. The above 
form, however, seems to require the fewest algebraic manipulations in the fol- 
lowing analysis • 

r Jpon use of the Green* s Function of the wave operator in unbounded space, 
Eq . (13) gives 


7 



4w*(x, t )=-/—(*+ — M # . . ) 6(f ) 5(q )dydx - — / — M * 5(f ) 5(g) dy dx 

~ •' r n c n t(y) ~ c3t J rn ** 

_ V-l I $ n 5( f ) 5 (q ) d£ dr (14) 

where j = t)/3t. 

If the spatial variables y are transformed to Lagrangian variables n, and 
if the transformation t-k} is employed and the delta functions integrated out, 
then Ea . (14) can immediately be written in the form of Eg. (5.3) of Reference 

[10]. Such an approach leaves the divergence and the time derivatinve oper- 
ators outside the integrals, however, which renders the resulting formula un- 
suitable for numerical application. In the current paper Eg. (14) will be 
manipulated differently to lead, as will be seen, to a more practically useful 
resu 1 1 . 

Before interpreting the integrals in Ea. (14), the following relation is 
used to convert the sDace derivatives in the last integral to a time 
derivative : 


3 .5(g), _ ! 3 r i<5(g) r i 5(g ) 

3xi r J ‘ c 3t l r J r 2 


(15) 


where r = (x-y)/r is the unit radiation vector. Note that in this relation the 
source variables (y,x) are held fixed. Eauation (14) can now be written as 


4tt$( 3 


x,t) = -/ ^-( $ n + — M n # x ( j ) 6(f ) 5(q )d£ dx + j - — $ cos 0 6(f) 6(g )d£ dx 


+ — “r / --(cos 9 - M ) $ 6(f )6(g)dy dx 
c dt 4 r n ~ 


(16) 
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There are two types of inteqrals in Eg. (16), which will be denoted 

1 1 = / 0 1 6 ( f ) 5(q )dy di (17) 

1 2 = 1 / Q 2 6(f)5(q)dydT (18) 

The technique of evaluation of these integrals will be given below. 

Although in many applications of the Kirchhoff formula in acoustics the 
surface S is taken as riqid, for completeness S is here assumed to be deform- 
able. Let S be described by the surface coordinates (u 1 , u 2 ) and assume 
u 3 — f • This mapping is only required in the vicinity of f=0. Let the mapping 
(u 1 , u 2 , u 2 ) > y be a differentiable function of time x* Denote the determi- 
nant of the coefficients of the first fundamental form on S by g(2) 3=1 ^11^22 ” 
q 12 2 , where the g— (i,j = 1,2) are the metric tensor components. Then g^ 2 ) i- s 
a function of time t in addition to (u\u 2 ). For example, consider a sphere of 
radius R(x) whose center is moving at speed U(x). Let (u^u 2 ) be the spherical 

surface coordinates on this sphere as shown in Fig. 1. Then if at x=0 the 

center of the sphere is at the origin of the £-frame, it follows that 

= y^(x) +[u 3 +R( t) ]cos u 1 sin u 2 

y 2 = + tu 3 +R( x) Isin u 1 sin u 2 (19) 

y^ = y^(x) + [ u 3 +R( x ) ] cos u 2 

where 

x 

y°(x) = / U(t') d x 1 (20) 

0 ~ 
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Also, it is seen that on S 


/q~ ( " 2) = R 2 d) sin u 2 


(21 ) 


To evaluate 1 1 , the followinq successive transformations are carried out: 

£ + (ui f u2, u 3) and then t + g. The Jacobians of transformations are 1 and 
1/(1-M r ), respectively, where = v *r/c . Here v = 3^/3 t holding (u^u^u 2 ) 
fixed. After the first transformation and integration with respect to f(i.e. 
u ^ ) , x j i s qi ven by 

t 

I = / / 0 (u 1 ,u 2 ,0,t) 6(g) /gTTT du 1 du 2 di (22) 

-oo d(S) [ 


1 2 

where D(S) is the domain of S in (u ,u^)-space. Note that even if S is deform- 
able, the limits of the inner integrals can be made time independent and are 
assumed so here (see Fiq. 1). The deformation of S appears in the variation of 
q ^ 2 ) with t* The second transformation t -► g and the subsequent integration 
with respect to g restrict the integrand to q = 0 and yield 


I 


1 


r 

I 

D(S) 


till , 

1-M 

r q=0 


du *du 2 


(23) 


The interpre tation of the inteqrand is as follows. If the relation g=0 is ex- 
plicitly written (remembering that u^=f=0 also in the integrand) 




t - t + x - y (u 1 ,u 2 ,0,t) /c = 0 , 


(24) 


then the solution in source time, t* , of this equation is the emission time 
which must be used in the integrand of Bg . (23). Note that since |v| < c by 
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assumption, Ba . (24) has a unique solution. Equation (23) will thus be written 

as 


I 


1 


/ 

D(S) 


0, ^ , -n 

— -] duMu 2 . 

i — M 


(25) 


Since D(S) is independent of t, I 2 follows immediately from Bg. (18) as 


1 , 3 ^ 2^(21 

c j at 


D(S) 


1-M 


du *du 2 = — 
I* c 


/ (■ 

D(S) 


___ 3__ Q 2 /g (2) .. 

M 3t * 1-M '' 

r r 


du *du 2 


(26) 


where the relation 3/ 3t = (1-M ) " 1 9/ 3 t | , obtained from Eq. (24), is used in 

/v 

the last expression. Note carefully that in taking 3/3 t the variables and 
are kept fixed. 

After using the results of Eqs. (25) and (26), Eg. (16) can be written as 


4ir#(x, t ) = -/ ( 

D ( S ) 


( t + c 1 M $ , ) /g. _ . 

n n T(y) (2) 


r(1-M r ) 


$/g cosQ 

du*du 2 + / [ ] du*du 2 

T D ( S ) r 2 ( 1 -M ) T 


‘■J 

D ( S ) 


( cos 9-M ) */g 

3 r n ( 2 ) ^ » 

_ — — y ~rk Ct \ J 1 . 


1-M 3t 


r ( 1 


-M ) 
r 


du *du 2 


(27) 


This result is equivalent to that of Morqans [5] and is referred to here as the 

qeneral Kirchhoff formula. The notation of the present paper, however, differs 

from Morqans’, who addressed the interior problem. 

There is a somewhat subtle point which should be made, however, about the 

result qiven in Reference [5]. The quantity [ /g ] du^u 2 is there written 

\ 2 ) T 

dS, but this is not the element of area of the physical surface S which is 
defined by fixinq x in f(y,T)=0. That element is /q ~~~ du^du 2 . The difference 
arises because as S deforms in space any element specified by (u ,u z ) has dif- 



ferent emission times t* for a fixed x and t. Corresponding to each t* the 
surface S has a different shape. Although the main result of Morgans is 
correct if interpreted as in the present paper, in some of his elementary 
examples there is an ambiguity in the meaning of dS . It should be noted that 
Morgans* result (27) is claimed to be erroneous by the authors of Reference 
[10]. However, they do not specifically point out the nature of the error. 

The modern derivation given above indicates that Morgans* formula is, in fact, 
correct, a conclusion which was made also by Munro [11]. 

Finally, it must be mentioned that the moving surface S in the above 
analysis can be unbounded provided that it divides the space into interior and 
exterior parts. An open bounded or unbounded surface can be considered as the 
limiting case of a closed surface of infinitely small thickness to which the 
Kirchhoff formula can be applied. 

If S is stationary, then M n and M r are zero and q^ 2 ) is time independent. 
Also x* is simply t-r/c. Under these conditions, Eq. (27), reduces to the 
Kirchhoff formula for stationary surfaces, Eq. (10). 


ALTERNATE FORM FOR APPLICATIONS 

Because Eq. (27) contains $ , v and also a time derivative of a compli- 

t ( y ) 

cated function in the last inteqrand, it is still not in a form which is 
readily applied for numerical analysis. In order to obtain a result which is 
practically useful it is essential that the time derivative in the integrand of 
the last term be evaluated analytically and the entire result written in terms 
of qeometric and kinematic quantities which are available or can be readily 
measured. To this end the following relations will be used: 


1 *L _ L L L 

c 9 T C 3 T ~ 


y(u ! ,u 2 ,0 , t) 


-M 


r 


(28-a) 


1 2 
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1 I'll = 1 M • r + J-(m 2 -m 2 ) 
c 3x c - - r r 


( 28-b) 


1 8cos0 1 • 1 / u u v 

n • r + — ( M cos 9-M ) 

c 8 t c ~ ~ r r n 


(28 -c) 


1 = — (A • n + M • n ) 

c9t c ~ ~ ~ ~ 


( 28-d ) 


i y 

where the dot over M and n denote source time derivative keeping (u ,u *) 
fixed. In addition, the following definitions are introduced: 


A = A • r 
r ~ ~ 


(29-a) 


n = n • r 
r ~ ~ 


(29-b) 


A = A • n 
n ~ ~ 


(29-c) 


n = n • M 


^ 'V 


( 29-d ) 


1 9 /<? ( 2 ) 


/q 


3t 


= a 


( 29-e ) 


( 2 ) 


Note that in the last inteqral of Eq. (27), a new time rate of chanqe of * will 

appear when (u 1 ,!) 2 ) is kept fixed. This is denoted *. The relation between 

4> , . and $ is obtained as follows: 

x(y ) 

,9 . . r 3* 3y i 3*, 

* = {37 *[y(u,T),T]} U3=0 = {37- 77- + a7> U3=0 


= fv • V$ + $ , , } = c M $ + cM • 

~ x( y ) u^=0 n n ~t 


V + *T(y) 


1 3 


( 30 ) 


where u = (u ] ,u 2 ,u^) f v = 9^/3r| u and is the tangential component of the 
Mach number vector M = v/c on S. Also stands for the surface gradient 
operator . 

Equation (27) can now be written as 


^1 ^ ( o ) /q / \ 

4 tt $( x , t ) = / [ — r- — - — r] du*du 2 + / [ H ^du^u 2 

D(S) r(1 “V T ' - 


J l J c 

D ( S ) r 2 (1-M ) 


(31 ) 


where the expressions E 1 and E 2 are 


*] * ™" 1 • • • 

E = (M 2 -1)$ + MM*V<l>-c M $ + — [ ( n - M - n )$ 

1 n n n~t 2 n 1 -M r n M 


* • r ' 1 

+ (cos9-M ) $ + (cos^-M ) 4>a] + — 1 

n n ( 1 -M )2 r 


[M ( COS 0-M )$] ; 

" n 


( 32-a ) 


E = cos 9 + — ( 2M cos 0-M M -M ] + ((cos 0-M )(M 2 -M 2 )] ( 32— b) 

2 1_M r r 0 r n (1-M ) 2 n r 

r 


If the surface S is assumed rigid, then 9(2) inde Pendent °f T and 


n = (jl i x n 


(33) 


where u> i s the angular velocity of S. One can now write dS = /g ~~~ du^du^ for 
the element of area of S. Equation (31) can thus be written for a rigid sur- 
face as 


E 1 

47r*(x,t) = / [--(fliTTh* dS + / (— J 


dS 


(34) 


S r 2 ( 1 -M ) 


where E 1 and E 2 are again given by Eg . (32) using Eq • (33) where needed. The 
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general Kirchhoff formula of Eqs. (31) and (34) is now in a form which is suit- 
able for practical applications. These results appear not to have been pub- 
lished elsewhere. Once more it can be checked that the Kirchhoff formula for 
stationary surfaces, Eq . (10), can be obtained from these equations. Again it 
can be seen that both Eqs. (31) and (34) remain valid for a piecewise smooth 
surface S. 


REMARKS ON EXTENSION TO SUPERSONIC MOTION 

A more general Kirchhoff formula can be derived which is valid for deform- 
able, piecewise smooth surfaces in arbitrary motion (subsonic or supersonic). 
The derivation of this result is very lengthy and complicated, and it will not 
be written out here. It relies on the solution of the following wave equation 
given by the authors [12]: 


□ 2 * = *{£l vf l <5(f ) ] . 


(35) 


In this equation = ( V, c 1 9/ 3t ) and g is a 4-vector. Note that |vf|=1 in the 
present paper. Equation (13) can be written as 


□ 2 * = -(* + — M $ )5(f ) - ? -[Q 6(f)] 

L ~ J n c n t(x) 4 ~ 


(36) 


where Q = <t>(n,M ). Now the Green 1 s function of the wave equation can be applied 
~ ~ n 

to find 4> utilizing the result of reference [12] for the second term of Eq. (36). 
The resulting Kirchhoff formula is more general than the formula presented here. 
However, the various forms of the general Kirchhoff formula of the current paper 


appear to be adequate for present applications in acoustics 



APPLICATIONS IN ACOUSTICS 


As mentioned in the introduction , application of the Kirchhoff formula has 
been proposed for rotor noise prediction. Attempts have also been made by 
Forsyth and Korkan [13] to use this formula for high speed propeller noise 
calculation. These authors have proposed using the Kirchhoff formula for a sta- 
tionary surface which, as shown above, is very different from the general 
Kirchhoff formula. The surface S in the work of Forsyth and Korkan is taken as a 
finite cylinder, fixed with respect to the propeller, with its axis along that of 
the propeller. This surface is therefore in motion with respect to the unbounded 
acoustic medium, and the general Kirchhoff formula must be applied to S. 

One of the interesting applications of the general Kirchhoff formula may be 
in the prediction of the noise of high bypass ratio turbofans. These propulsors 
are currently favored in aircraft designs because of their high efficiency. By 
using a cylindrical surface S with its axis along the fan axis and its ends at 
the fan inlet and exhaust, Fig. 2, the main result of the present paper can be 
applied to predict the radiation pattern. There are several sophisticated turbo- 
machinery aerodynamics codes which can supply the needed input data to the 
Kirchhof f formula. 

Lyrintzis and George have developed a specialized Kirchhoff formula to cal- 
culate the far field noise from a single helicopter blade interacting with a 
vortex [14]. Their formula applies to uniform motion of a rigid surface and is 
derived by using Green's theorem. It is equivalent, in their case, to the 
general Kirchhoff formula derived here, Bg. (34). 
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Figure 2: Illustrating a surface S appropriate for 

turbofan noise analysis. 
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